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Abstract 

Given relatively prime positive integers ai,..., a^, the Frobenius number is the 
largest integer that cannot be written as a nonnegative integer combination of the 
Oj. We examine the parametric version of this problem: given a* = Ui (t) as functions 
of t, compute the Frobenius number as a function of t. A function / ; Z_|_ —)■ Z is 
a quasi-polynomial if there exists a period m and polynomials /o,..., fm-i such 
that f{t) = ft mod m(d for all t. We conjecture that, if the ai{t) are polynomials (or 
quasi-polynomials) in t, then the Frobenius number agrees with a quasi-polynomial, 
for sufficiently large t. We prove this in the case where the ai{t) are linear functions, 
and also prove it in the case where n (the number of generators) is at most 3. 


1 Introduction 

Given positive integers a*, 1 ^ i ^ n, let 

Pi £ ^>0 

be the semigroup generated by the Oj. If the a* are relatively prime, dehne the Frobenius 
number F{ai,..., a„) to be the largest integer not in (oi,..., a„). The Frobenius prob¬ 
lem of determining F(ai,..., a^) has a long history — Sylvester proved [8] in 1884 that 
F{a,b) = ab — a — b, and see Ramirez Alfonsin’s book [5] for many subsequent results. 

It will be convenient to also dehne F(ai,..., a^) in the case where the Oj are not 
relatively prime, so that there is no largest integer not in the semigroup. A reasonable 
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definition seems to be the largest integer in the group Zoi + ■ ■ ■ + that is not in the 
semigroup (oi,..., a„). That is, if d is the greatest common divisor of oi,..., On, then 
F{ai,..., Gn) = dF {^,..., ^). Sylvester’s identity then becomes F{a, h) = lcm(a, h) — 
a — b. 

The parametric Frobenius problem is, given functions Oi : Z_|_ —)■ Z_|_ to determine 
F(ai(t),..., Gnit)) as a function of t. For example, using Sylvester’s identity, 

^ ^ ^ ^ — t — (t + 2) if t is even. 

We see that F(t,t + 2) is a quasi-polynomial] a function / : Z+ —)■ Z is a quasi¬ 
polynomial if there exist an m G Z+ and polynomials /o,...,/m-i G Q[f] such that 
fit) = ft mod mit) for all t G Z+. Here m is a period of / and the ft are components of /. 
We will assume that all of our functions are integer-valued, but note that the F’(t, t -|- 2) 
example shows that we may need polynomials with rational coefficients. 

In general, our functions may misbehave for small t. We say that a property is even¬ 
tually true if it is true for all sufficiently large t. We say that fit) G EQP (short for 
eventual quasi-polynomial) if fit) eventually agrees with a quasi-polynomial. We say 
that the degree of / G EQP is the maximum degree of its components. 

We conjecture that, if aiit) G EQP, then F(ai(t),... ,a„(t)) G EQP. We prove this 
for the special case where the Giit) are linear functions and for the special case where 
n ^ 3. 

Conjecture 1.1. Suppose that aft) G EQP, 1 ^ i ^ n, are eventually positive. Then 
F{aiit),...,anit)) gEQP. 

Theorem 1.2. Suppose that Giit) G EQP, 1 ^ i ^ n, are eventually positive and have 
degree at most 1. Then F(ai(t),..., a„(t)) G EQP, with degree at most 2. 

Theorem 1.3. Suppose that aft) G EQP, 1 ^ i ^ n, are eventually positive, with n ^ 3. 
Then F(ai(t),..., Gnit)) G EQP. 

These results are examples of so-called “unreasonable” appearances of quasi-polyno¬ 
mials, as discussed by Woods [10]. “Reasonable” appearances trace back to Ehrhart’s 
classical result [3] that, if P C M"- is a polytope with rational vertices, then fit) = 
ifitP n IT) is a quasi-polynomial. Note that if P is dehned with linear inequalities 
bj ■ X ^ Cj, then tP is dehned with linear inequalities bj ■ x ^ cR. As t changes, these 
linear inequalities move, but their normal vectors (bj) remain the same. Indeed, Woods 
proved [11] that any example over the integers dehned with linear inequalities, boolean 
operations (and, or, not), and quantihers (V, 3) has this quasi-polynomial behavior. This 
is true even if there is more than one parameter; for example, 

#{ix,y) eZlQ\2x^t,3y^s} = (^ ^ + 1 ^ + 1 ^ 

is a quasi-polynomial of period 2 in t and period 3 in s. 
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If we look at the parametric Frobenius problem, however, we see that it does not fit 
this pattern. For example, if we want to ask whether u G (t, f + l,t + 2), we are asking 
whether the polytope 


{(x, y, z) e M|o \ u = tx + {t + l)y + {t + 2)z] 

contains any integer points. For a hxed t, this is a 2-dimensional triangle in as t 
changes, this triangle “twists” (the normal vector changes). Examples such as this were 
categorized in [10] as “unreasonable”, though they are conjectured to still lead to eventual 
quasi-polynomial behavior. 

This paper adds a third example to the list of recent results demonstrating this phe¬ 
nomenon; previously Chen, Li, and Sam [2] showed that the number of integer points in a 
polytope whose vertices are rational functions of t is in EQP, and Calegari and Walker [1] 
showed that the vertices of the integer hull of such a polytope have coordinates in EQP. 
A critical tool used in all of these results is that the division algorithm and the gcd of 
polynomials has quasi-polynomial behavior (cf. Lemma 3.1); for example, the Euclidean 
algorithm yields 


gcd(2t -I- 1, 5t -I- 6) = gcd(t -I- 4, 2t -I- 1) = gcd(7, t -|- 4) 


7 if t = 3 mod 7, 
1 otherwise. 


Note that unlike in the “reasonable” case, these results only hold with one parameter 
variable. For example F{s, t) = lcm(s, t) — s — t is not a quasi-polynomial in s and t. 

The original inspiration for this paper comes from a conjecture of Wagon [9] (see [4, 
Section 17]) that, for any hxed M and residue class j of t mod M^, there exist (usually 
positive) integers CM,j and dM,j such that eventually 

F(t, t + ,t + 2 '^,... ,t + M^) = {t^ -\- CM,jt^ — dM,j- 

Using the proof of Theorem 1.2, we can prove that there is indeed quasi-polynomial 
behavior with period M^. In fact, such a result holds in greater generality: 


Corollary 1.4. Suppose bi < b 2 < <■■■< bn are integers. Then F(t + bi,t + b 2 ,... ,t + 

bn) G EQP with period bn — bi. 


When the bi form an arithmetic sequence, a precise quasi-polynomial formula was 
previously given by Roberts [6]: 


F(t ,t + d.,... ,t + sd) 


( 

t-2 


s 



t + {d-l){t-l) - 1. 


In Section 2, we work through the example 


E(t, t -|- 1, t -|- 2) 


t-2 



t-l = 



if t even, 
— 1 if t odd. 
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which will give a flavor of our general proof. In Section 3, we prove Theorem 1.2. In 
Section 4, we prove the various lemmas needed. In Section 5, we prove Corollary 1.4. 
In Section 6, we prove Theorem 1.3, using Rpdseth’s algorithm [7] for the 3 generator 
Frobenius problem. 

The proof of Theorem 1.2 relies on the fact that semigroups with only two generators 
are much easier to deal with. Indeed, we have the following definition and lemma: 

Definition 1.5. Let a, 6 G be relatively prime, and let c G Z. The canonical form for 
c is given hy c = pa + qh with p,q E Z and 0 ^ p < b. 

Lemma 1.6. Let a,b E be relatively prime, and let c E Z. 

1. The canonical form for c exists and is unique. In particular, if c = p'a + q'b is 
any form with p',q' E Z, and if k and r are the quotient and remainder when p' is 
divided by b, then the canonical form for c is ra+ {q' + ka)b. 

2. If c = pa + qb is in canonical form, c G (a, b) if and only if q ^ 0. 

2 An Example 

In this example, we compute F(t, t + l,t + 2). Let a = t, b = t + 1, and c = t + 2 be our 
three generators, let S = (a, b, c), and let T = (a, c). Notice that 26 = a + c. This implies 
that S = T U (6 + T), as follows: if pa + qb + re is a representation of an element of S, 
with p,q,r E Z^q and q ^ 2, then (p + l)a + {q — 2)6 + (r + l)c is also a representation, 
so we may assume without loss of generality that g is 0 or 1 (cf. Lemma 3.4). Next, we 
run the extended Euclidean algorithm on the integers a and c, and get 

gcd(a, c) = gcd(f, t + 2) = gcd(2, t) 

(cf. Lemma 3.1). The next division step in the algorithm depends on the parity of t, and 
rather than end up with the messy [t/2j, we divide into two cases. 

The case t is odd: Let t = 2s + 1, so that a = 2s + 1, 6 = 2s + 2, and c = 2s + 3. 
Now 

gcd(a, c) = gcd(2,2s + 1) = gcd(l, 2) = 1, 

and the extended Euclidean algorithm yields 1 = (s + l)a — sc. Let u E Z, and we are 
wondering whether u ^ S, that is, u ^ T and u ^ b + T. Suppose that 

u = pa + qc 

is the canonical form for u (see Definition 1.5). Lemma 1.6(2) tells us that u eT ii and 
only if g ^ 0. Since we are looking for u ^ T, we may assume from here on out that 
g < 0. To characterize when u E b + T (that is, u — b E T) we must And the canonical 
form for u — b. 
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First we compute the canonical form for b. Multiplying the equation 1 = (s + l)a — sc 
by 6 = 2s + 2 yields some form for b: 

b = (2s + 2)(s + l)a — (2s + 2)sc. 

By Lemma 1.6(1), we may hnd the canonical form by dividing (2s + 2)(s + l) by c = 2s + 3: 
the quotient is s with remainder s + 2, giving the canonical form for b as 

6 = (s + 2)a — sc. 

Note that we got a little lucky here: our remainder of s + 2 is clearly less than c = 2s + 3; 
if our remainder had instead been s + 7, say, then we would only have the canonical form 
for sufficiently large s. 

Now we have that some form for u — b is 

u — b = {p — s — 2)a + {q + s)c. 

Is this the canonical form? There are two cases: 

If p ^ s + 2, then O^p — s — 2^p<c, and this is in canonical form. Therefore, to 
have u — b E T, we must have g + s ^ 0 (again using Lemma 1.6(2)), that is, q ^ —s. 

If p < s + 2, then the canonical form is 

u — b = (p — s — 2 + c)a + (g + s — a)c = (p + s + l)a + {q — s — l)c 

(to check that this is canonical, note that p ^ 0 and s + 2 < c imply p — s — 2 + c^O, 
and p < s + 2 implies p — s — 2 + c<c). In this case, to have u — b E T, we must have 
g ^ s +1. This implies that g ^ 0, but we have assumed g < 0 (so that u ^ T). Therefore 
this case never has u — b E T. 

To summarize across both cases, if u = pa + gc, with 0^p<c = 2s + 3in the 
canonical form, then u ^ S' if and only if 

not(g ^ 0) and not((p, g) ^ (s + 2, —s)) 

(cf. Lemma 3.5). The set of such (p, g) has a “stairstep” shape, as can be seen in Figure 

1 for t = 5. In particular, the set of such (p, g) can be rewritten as p ^ 0 and 

(p ^ s + 1 and g ^ —l) or (p ^ 2s + 2 and g ^ —s — l) 

(cf. proof of Lemma 3.6, d = 1 case). The two “corners” in this picture, where both 

inequalities of one of these two conjunctions are tight, give our candidates for the largest 
u ^ S: it must be either 

(s + l)a — Ic = 2s^ + s — 2 or (2s + 2)a + (—s — l)c = 2s^ + s — 1. 

The latter is always larger (in general, one might only be eventually larger than the other), 
and so we have ^ 

F{t, t + l,t + 2) = 2s^ + s — 1 = — — - — 1, 
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Figure 1: The t = 5 case of the example. We have T = (5, 7) and S = (5, 6, 7) = TU(6+T). 
Points {p, q) G with 0 ^ p < 7 are the canonical forms for u = 5p + 7g G Z. Then 
g ^ 0 corresponds to n G T, (p, q) ^ (4, —2) corresponds to n G 6 + T, and positive u 
such that u ^ S are labelled beside their corresponding (p, g). The “corners” u = 8 and 
u = 9 are candidates for the Frobenius number, and so F{5, 6, 7) = max{8, 9} = 9. 


in this case where t = 2s + 1 is odd. 


The case t is even: Let t = 2s, so that a = 2s, 6 = 2s + 1, and c = 2s + 2. To 
have M ^ S' we must have both u ^ T = {a, c) and u ^ b + T. Note that every element 
in T is even and every element of 6 + T is odd. Therefore the largest even integer not 
in S is the largest even integer not in T, which since T has only 2 generators is simply 
lcm(a, c) — a — c = 2s^ — 2s — 2 (we’re getting lucky here that we didn’t have to do an 
analysis similar to the t is odd case). Similarly the largest odd integer not in S is the 
largest odd integer not in 6+T, which is 6+(2s^ —2s —2) = 2s^ —1. The largest integer not 
in S is then the maximnm of these two candidates, so F{t,t + l,t + 2) = 2s^ — 1 = b- — 1 
(cf. proof of Lemma 3.6, d > 1 case), in this case where t is even. 


Combining the even and odd case, we have proved that 


F(t, t + 1, t + 2) 



if t even, 
— 1 if t odd. 


3 Proof of Theorem 1.2 

This hrst lemma gives us the basic tools we need: 

Lemma 3.1. Given f,g G EQP, 

1. There exists an m E Z+ such that, for each 0 ^ j < m, f{ms + j) eventually agrees 
with a polynomial in Z[s]. 
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2. If deg{g) > 0, there exists q,r E EQP sueh that fit) = q(t)g(t) +r(t) and deg(r) < 
deg( 5 f). Note that this is the analogue of the traditional division algorithm over<Q\f\. 

3. If git) is eventually positive, there exists q,r E EQP sueh that fit) = qit)git) +r(t) 
and eventually 0 ^ r(t) < git). Furthermore, deg(r) ^ deg( 5 f). Note that this is the 
analogue of the traditional division algorithm over Z. 

4- There exists p,q,d E EQP such that gcd (^f{t),g{t)) = d{t) and d{t) = p{t)f{t) + 
q{'t)g{t). 

5. We have max (^f{t),g{t)) E EQP. 

These are proved in Section 4 of [2] (and also in [1]), so in Section 4 we merely give 
an outline of the important steps. 

Remark 3.2. Lemma 3.1(1) will allow us to often simply say “without loss of generality, 
/ G Z[f]”: we may analyze f{ms + j) for each j, recognize that statements may be false 
for small s, and then convert back to t using s = {t — j)/m. 

Using this remark, we may assume that the ai{t) are polynomials in Z[t]. Since they 
are of degree at most 1 and eventually positive, we have that 


Qiit) = ait + Pi, 

with either ctj G Z+, G Z or o;* = 0, Pi E TL^^. Furthermore, without loss of generality, 
we may assume that they are ordered so that 

Piaj ^ Pjai 


for i ^ j, that is (for ai ^ 0), 

A < ^ < Pn 

ai a2 a^i 

We hrst consider a degenerate case, where Pia^ = P^ai, for all i,fj. 

Lemma 3.3. If Piaj = Pjai, for all i,j, then there are ao, Pq E Z^q such that, for all i, 
ai = 7jQ;o ond Pi = 'jiPo (for some y* G Z+j. Therefore 

F{ai{t),Unit)) = F( 7 i, ..., 7n) ■ {aot + Po) 

is a polynomial of degree at most 1. 

Therefore, we can assume that Pian < Pndi- In particular, the polynomials ai and 
On do not share a linear factor, and so if d{t) = gcd [apt) , an{t)), then d{t) E EQP is 
of degree 0, that is, dit) is eventually a periodic function. Let Sit) = {apt),... ,a„(t)), 
and let T{t) = (apt), an{t)). T{t) is a semigroup with only two generators, so it is much 
easier to analyze. The following lemma allows us to cover S{t) with a hnite number of 
translated copies of T{t). 
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Lemma 3.4. For ai(t) = aF + f^i as described above, there exists a finite set H of integer¬ 
valued polynomials of degree at most 1 such that 

(ai(t),..., anit)) = IJ (^h{t) + (ai(t), a„(t)) j 

h&H 

(for t sufficiently large so that afit) > Offii). Furthermore, 0 G and the other h & H 
are eventually positive. 

Assume for the moment that afif) and an{t) are relatively prime. By Lemma 1.6(2), 
the set of all integers not in T will be in bijection to the set 

[{p,q) ^p < anit), g < 0} , 

under the bijection {p,q) e-)■ pafit) + qanit). By Lemma 3.4, integers not in S must also 
not be in h{f) + T, for all h & H. The following lemma gives an easy way to check this. 

Lemma 3.5. Let f,g,h G EQP, with fit) and g(t) relatively prime for all t and with 
fit), git), hit) eventually positive. There exists r,s E EQP such that, given p,q E Z with 
0 ^ p < git) and q < 0, 

pf{t) + qg{t) E hit) + {f{t),g{t)) if and only if {p, q) ^ (r(t), s{t)) 

component-wise. Furthermore, deg(r) ^ deg(g) and deg(s) ^ max{deg(/), deg(/i) — 
deg(g)}, 

In our case, this gives r, s of degree at most 1. If afit) and Onit) are relatively prime, all 
that remains is to follow the implications of Lemmas 3.4 and Lemma 3.5, which will give 
us a staircase-shaped set such as in Figure 1. If aiit) and a„(t) are not relatively prime, we 
must reduce to the case where they are. Both of these are accomplished in the following 
lemma, which (combined with Lemma 3.4 and the fact that d{t) = gcd (/(t),g(t)) is of 
degree 0) proves our theorem. 

Lemma 3.6. Let fit),git) E EQP, and let d{t) = gcd (/(t),g(t)) be of degree 0, let 
H C EQP be a finite set, and let 

S{t) = IJ {h{t) + {f{t),g{t))). 

h&H 

The function Fit) giving the largest integer not in Sit) is in EQP, with degree at most 
max{deg(/) -h deg( 5 (), deg(h)}. 

In our case, this gives that F is of degree at most 2. 



4 Proofs of Lemmas 


Proof of Lemma 1.6. Part 1 is a standard result, using the extended Euclidean algorithm 
and the fact that if c = p'a + q'b is one solution, then all solutions are given by c = 
(p' — kb)a + (g' + ka)b for k e Z. 

For part 2, we have c = pa + qb with 0 ^ p < b. The reverse implication is immediate: 
if g ^ 0 , then we have written c = pa + qb with p,q E Z^q, proving that c G {a,b). 
Conversely, suppose c G (a, b) so that c = p'a + q'b with p', q' G Z^q. By part 1, if k and 
r are the quotient and remainder when p' is divided by b, then 

c = ra + {q + ka)b 

is the canonical form for c, with q = q' + ka ^ Q. □ 

Proof of Lemma 3.1. As these are proved in Section 4 of [2], we simply give an outline 
here. 

For part 1, this is mostly obvious: simply take component polynomials of the quasi¬ 
polynomial. The main subtlety, as seen in the example from Section 2, is that integer¬ 
valued polynomials may have non-integral coefficients. In this case, let m be the least 
common multiple of the denominators of the coefficients. Examining f{ms + i), we see 
that all coefficients of must be integral, except possibly the constant coefficient; since 
the function is integer-valued, the constant coefficient must also be integral. 

For part 2, simply perform polynomial division. The main subtlety is the following: 
Suppose, for example, that f(t) = 1“^ + 3t and g(t) = 2t + 1. Then the leading coefficient 
of g does not divide the leading coefficient of /, and the traditional polynomial division 
algorithm would produce quotients that are not integer-valued. Instead, we look sepa¬ 
rately at each residue class of t modulo the leading coefficient of g; for example, if t is 
odd, then t = 2s -|- 1 for some s G Z^q, so substituting gives /(2s -|- 1) = 4s^ -|- 10s -|- 3 
and g(2s -|- 1) = 4s -|- 3, and now the leading term does divide evenly. 

For part 3, hrst perform the polynomial division algorithm for part 2. For example, 
suppose fit) = 2t — 3 and g(t) = t, yielding f = 2g -\ —3. For part 3, however, we want 
the integer division algorithm: fit) = lg(t) -|- (t — 3), and the remainder t — 3 is between 0 
and g as long as t ^ 3. In other words, if we have found / = q'g + r' with deg(r') < deg(g), 
but we eventually have r'(t) < 0, then we should use quotient g = g' — 1 and remainder 
r = g + r' instead, as eventually 0 ^ git) + r'(t) < git). 

For part 4, run the extended Euclidean algorithm, repeatedly performing the division 
algorithm of part 2 . 

For part 5, note that given two polynomials, one eventually dominates the other. □ 

Proof of Lemma 3.3. For each / let y* = gcd(Q;i,/3i) G Z+, let and let /3' = 

A/bi) with a' and 13'^ relatively prime. Dividing /3iaj = /3iai by 7 ^ 7 ^ yields Aa/ = 

Since a'- divides and is relatively prime to /S', we have that a'j divides a'. Similarly, 
a' divides a'-, and so a' = a'-. Similarly, A = Pj- Taking ao to be the common a' (equal 
for all i) and A to be the common A; the proof follows. □ 
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Proof of Lemma S.f. Let r = /9 iq;„. We are given that r > 0 and that /3jai—/3iaj ^ 

0 for all i ^ j. Let 


H — \^X20,2 + ■ ■ ■ + Xn-lO-n-l \ \ ^ ^>0) Xi < . 

We will prove that H has the required properties; certainly every element of H is of degree 
at most 1. We must show that, if m G (ai(t),..., an(t)), then there exists a A G such 
that u = XiOiit) + ■ ■ ■ + Xnanit) and A* < r for 2 ^ i ^ n — 1 . 

By dehnition of the semigroup, there exists a A G Z”g such that u = Aiai(t) + ■ • • + 
Xnanit). Suppose that Xi ^ r for some i with 2 ^ i ^ n — 1. Let p = — Attn and 

q = I3iai — I3iai and observe that p, g ^ 0. Then 

pai(t) + qan{t) = {Pnai - Attn)(ttit + A) + (Atti - Atti)(tt„f + A) 

= I3naiait - AttnA + Atti A - l^idiant 
= (Attl - Attn)(tt*t + A) 

= rai{t). 

Now dehne X' by A'^ = Ai +p. A' = Xi — r, X'^ = A„ + g, and A'- = Xj for all other j. Then 

u = X\ai{f) H-h AAn(t) 

is a new representation of u. Repeating this process eventually yields a representation 
with Xi < r for 2 ^ i ^ n — 1. □ 

Proof of Lemma 3.5. Assume 0 ^ p < g{t) and g < 0. We need conditions on p,q for 
which pf{t) + qg{t) — h{t) G {f{t),g{t)). By Lemma 3.1(4) and the fact that f{t) and 
g(t) are relatively prime, we may write 1 = r'(t)f(t) + s'(t)g(t), where P,s' G EQP. 
Multiplying this equation by h(t) yields 

h{t) = {h{ty{t))f{t) + {h{t)s'{t))g{t). 

If we let k(t) and r(t) be the quotient and remainder when h(t)r'(t) is divided by g(t), 
using Lemma 3.1(3), and let s(t) = h(t)s'(t) + k(t)f(t), then we have 

hy = r{t)f{t) + s{t)g{t), 

with 0 ^ r(t) < g(t). In particular, this gives us that deg(r) ^ deg(g) and therefore 
deg(s) ^ max{deg(/), deg(/i) — deg(g)}. We have 

Pfit) + qg{t) - h{t) = (p - r{t))f{t) + (g - s{t))g{t). (1) 

The case p ^ 'r(t): Then 0 ^ p — r{t) < g{t), and (1) is in canonical form. Then 
P/(A + qgy ~ hit) ^ {/(A)5'(^)) if only if g — s{t) ^ 0, by Lemma 1.6(2). 

The case p < r{t): Then the canonical form for pf{t) + qg{t) — h{t) will be 

(p - r{t) + g{t))f{t) + (g - s{t) - f{t))g{t), ( 2 ) 
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since 0 ^ p < r(t) and r(t) < g(t) imply that 

0 < p - r{t) + g{t) < g{t). 

Since h(t), r(t), f(t), and g(t) are eventually positive and r(t) < g(t), we eventually have 

0 < h{t) = r{t)f{t) + s{t)g{t) < g{t){f{t) + s{t)), 

and so eventually f(t) + s(t) > 0. Since we are assuming that q < 0, eventually q — s(t) — 
fit) < 0. Therefore the canonical form (2) shows that pfit) + qgit) — hit) ^ (fit), git)), 
by Lemma 1.6(2). 

Combining the two cases, we see that pfit) + qgit) G hit) + (fit), git)) exactly when 
p ^ r(t) and q ^ sit), as desired. □ 

Proof of Lemma 3.6. Since dit) G EQP is of degree zero, it is eventually a periodic func¬ 
tion. By Remark 3.2, we may focus on a component of dif) and assume that dif) = d is 
a constant. 

The case d = 1: We assume, without loss of generality, that 0 G iL and the other 
functions in H are eventually positive: if not, let h^it) be the eventually minimal polyno¬ 
mial in H, and hnd the largest integer Fft) not in 

S'it) = IJ ({hit) - hoit)) + {fit), git)) 

h&H h 

then Fit) = Fft) + hoit). 

We wish to describe the set U it) of (p, q) corresponding to canonical forms u = pfit) + 
qgit) such that u ^ S. Canonical implies that 0 ^ p < git), and u ^ iO + {fit), git))) 
implies that g < 0, by Lemma 1.6(2). Each h G iL \ {0} gives the condition “p < 
Thit) OT q < Shit)”, by Lemma 3.5. Then the set f/(t) has a “stairstep” shape as in 
Figure 1. To be precise, order the {rhit), Shit)) such that ri(t) ^ r 2 (t) ^ ^ 'rm(t) 

(eventually). Notice that if rj(t) ^ rjit), then we may assume without loss of generality 
that Siit) > Sjit) (eventually), or else the condition “p < Vjit) or q < Sjit)” would be 
redundant. Further include (ro(t), So(t)) = (0,0) (corresponding to m ^ 0 -|- {fit), git))) 
and {r^n+iit), Sm+iit)) = {git), —fit)) (corresponding to (p, g) being canonical). Finally, 
for 0 ^ i ^ m, let {aiit),/3iit)) = (rj+i(t) — l,Si(t) — l). Then the set of (p, g) G f/(f) is 
exactly the set such that 0 ^ p and 

(p,g) ^ (ao(t),/3o(t)) or ■■■ or (p, g) ^ (^^(t),/3m(t)). 

Since if (p, g) ^ {aiit), fdiit)), then pfit) + qgit) ^ aiit)fit) + /di(t)g(t), our candidates 
for the largest integer not in Sit) are the aiit) fit) + (3iit)git), 0 ^ i ^ m. Then our hnal 
answer is 

max aiit) fit) -F A(t)g(t), 

which is in EQP, by Lemma 3.1(5). Furthermore, using that deg(Q;j) = deg(rj) ^ deg(g) 
and deg(/dj) = deg(si) ^ max{deg(/), deg(/i) — deg(g)}, by Lemma 3.5, we have that our 
hnal degree is at most max{deg(/) -|- deg(g), deg(/i)}. 


11 



The case d > 1: Let Hj{t) = {h & H \ h{t) = j mod d}. The remainder when a 
given h(t) is divided by d is a periodic function, by Lemma 3.1(3); applying Remark 3.2 
as necessary, we may assume that these remainders are constant, that is, that Hj = Hj{t) 
does not depend on t. Let Sj(t) = {u E S(t) \ u = j mod d}. Since d divides every element 
of {f{t),g{t)), 

Sj{t)= IJ (h{t) + {f{t),g{t)) 

heHj ^ 

Let Fj(t) be the largest integer not in 




U 

heHi 


hit) - j 
d 


+ 


M 9M\\ 

d ' d / J' 


By the case d = 1 proved above, T)(t) is in EQP. Since dFjit) + j is the maximum 
u E j + dZ such that u ^ Sjit), we get that F(f) = max^ (dFjit) + j). Lemma 3.1(5) then 
implies that F(f) is in EQP. □ 


5 Proof of Corollary 1.4 

Without loss of generality, we may assume that 6i = 0 (substituting s = t — hi, if 
necessary). For a hxed j E Z, look at all t = j mod bn, and let s be such that t = bnS + j. 
Let ttiis) = t + bi = bnS + j + bi. We must prove that F(ai(s),..., a„(s)) is a polynomial 
in s. We follow the proof of Theorem 1.2 and of the lemmas, looking for places where 
periodicity might be introduced. Note that the ai(s) are correctly ordered (as dehned 
before the statement of Lemma 3.3) so that Lemma 3.4 gives us 

*5'= U (d(s) + (ai(s),an(s))). 

h&H 

Following the proof of Lemma 3.4, note that each h G id is a nonnegative integer 
combination of 02 ( 5 ),..., an-i(<s). In particular, each h E H has the form kbnS + i for 
some integers k and i. 

Focus on a specihc h = kbn-s + i- Let 

d(s) = gcd (ai(s), a„(s)) = gcd(6„s + j, bnS + j + bn) = gcd(&„, bnS + j) = gcd(j, bn), 
which is a constant. Finding p and q with pj + qbn = d gives us 

d = ips-q + p)aiis) - (ps - g)a„(s). 

Assume for the moment that d = 1. We next examine the proof of Lemma 3.5. The hrst 
step is to take 

h = h ■ 1 = h ■ [ips — q + p)aiis) — (ps — g)a„(s)) 

= {ikbnS + i)ips - q + p))aiis) - ((/c6„s + d)(ps - g + p))a„(s). 
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and write it in canonical form. This involves taking the remainder when dividing 

{kbnS + i){ps - q + p) = kpbnS^ + {^p - kbnQ + kbnp)s + {-^q + ^p) 

by a„(s) = bnS + j + bn- The first step of polynomial long division works over the 
integers, because the leading coefficient, 6 ^, of an{s) divides into the leading coefficient, 
kpbn, of the dividend. This leaves a remainder that is linear in s. One more step of long 
division (dividing a linear in s function by a linear in s function and taking the linear in 
s remainder) will then lead to the hnal answer, without introducing extra periodicity. 

No other steps in the entire proof have a possibility of adding periodicity, so we are 
done, in the case d = 1. 

If d > 1, following the proof of Lemma 3.6 in the d > 1 case, let b'^ = bn/d, j' = j/d, 
and i' = [i/d\ be integers. Then the reduction to the d = 1 case gives that we must 
examine, 

), 

and we have reduced to the d = 1 case. 

6 Proof of Theorem 1.3 

The n = 1 case is trivial: F(ai(t)) = is the usual interpretation. 

The n = 2 case follows from Sylvester’s formula: 

F(ai(t), a 2 (t)) = 1 cm (ai(t), a 2 (t)) - ai{t) - a 2 (t) = Qi(^)Q 2 (t) _ 

gcd(ai(t),a 2 (t)) 

and the gcd can be computed using Lemma 3.1. 

For the n = 3 case, we show that the steps of Rpdseth’s algorithm [7] can be performed 
on eventual quasi-polynomials. We enumerate the steps of the algorithm as described for 
integers, and follow each step with a comment on how it works with elements of EQP. 

(1) If d{t) = gcd (ai(t), a 2 (t), a 3 (t)) ^ 1, then use that 

For elements of EQP, Lemma 3.1(4) allows us to compute this gcd. 

(2) We may assume that ai(t), 02 (t), 03 (t) are relatively prime where ai(t), 02 (t), a^it) G 
EQP. If e{t) = gcd (ai(t), a 2 (t)) 7 ^ 1, then use that 

F{ai{t),a 2 {t), as{t)) = e{t)F (e(t) - l). 


h{.s) 


Cil(s) (ln(s) 


d 


= {kb'^s + i') + {b'^s + f, b'^s + i' + b'^ 
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This can again be done using Lemma 3.1(4). 


( 3 ) We may assume that ai(t) and 02 (t) are relatively prime. Compute So(^) ^ EQP 
such that a 2 (t)so(t) = a^it) mod ai(t) and 0 ^ So(t) < ai(t). 

To do this for elements of EQP, Erst use Lemma 3.1(4) to find p,q & EQP such that 
1 = p{t)ai{t) + q{t)a 2 {t). Then 

a2(t){a3(t)q(t)) = a3(t){q(t)a2(t)) = a^it) ■ 1 mod ai(t). 

Now let So(t) be the remainder when a 3 (t)q(t) is divided by ai(t), using Lemma 3.1(3). 
This is the desired So(t). 

(4) If so(t) = 0, then a^it) is a multiple of ai(t), and the Frobenius problem reduces to 
the n = 2 case. 

If some of the components of So(t) G EQP are zero, we will have to split into cases, one 
for each component. 

(5) Compute qi,Si G EQP such that 

«i(f) = qiit)soit) - 0 < Si{t) < So{t). 

If deg(so) = 0, this is a slight variant of the usual integer division algorithm. If deg(so) > 
0 , use Lemma 3.1(2) to compute ai(t) = q(t)so(t) + r(t), with deg(r) < deg(so). If r(t) 
is eventually negative or zero, simply set qi(t) = q(t) and Si(t) = —r(t), and we have 
that eventually 0 ^ si(t) < so(t). If r(t) is eventually positive, set qi(t) = q(t) + 1 and 
Slit) = Soit) — r(t), which will eventually satisfy the bounds. 

( 6 ) Continue computing 

Soit) = q 2 it)siit) - S 2 (t), 0 ^ S 2 it) < Slit), 

Slit) = q3it)s2it) - Ssit), 0 ^ Ssit) < S 2 (t), 

Sm—2it) — qmit)Sm—lit) Smit), 0 ^ Smit) >Sm—l(t)) 

'^m—1 it) = qm+lit) it) + <Sm+l — 0 . 

This is simply repeating the process from Step 5. We must establish that it terminates. It 
suffices to show that, if deg(sj) > 0 , then deg(sj+ 2 ) < deg(sj) (and then once deg(sj) = 0 , 
the integers s* strictly decrease, so this will eventually terminate). Indeed, as described 
in Step 5, Lemma 3.1(2) gives us Si_i(t) = g(t)si(t) + r(t), with deg(r) < deg(si). If r(t) 
is eventually negative or zero, we have Sj+i(t) = —'r(t), and the degree has decreased. If 
r(t) is eventually positive, we have Sj+i(t) = Siit) —'r(t), which still have the same degree 
as Si- But in the next step, Si+i(t) = 1 ■ Sj(t) — r(t), and so Si+ 2 (t) = r(t) has lower 

degree. Of course, different components of the eventual quasi-polynomials may terminate 

at different m, and we must analyze each separately. 
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(7) Define s_i(t) = ai(t), -P-i(t) = 0, Po(t) = 1, and recursively Pi+i(t) = qi+i(t)Pi(t) — 
Pi-i{t). Then 

Srn+^ Sr^ So{t) g_i(t) _ 

Pm+l{t) Pm{t) ■" Po{t) P-l{t) 

Determine the unique index i such that 

■Sj+i(t) ^ asjt) ^ gj(t) 

Pi+i(t) ^ a2(t) Pi{t)' 

Restricting to a hxed residue class modulo the period of the quasi-polynomials, such an 
index i exists, since each Si{t)/Pi{t) is a rational function, so is eventually greater than, 
eventually less than, or eventually equal to ay,{t)/a- 2 {t). 

(8) Then 

F(ai(t),a2(t),a3(t)) = -ai{t) + a2{t) (sj(t) - 1) -h a^it) {Pi+i{t) - l) 

— min {a2{t)si+i{t),as{t)Pi{t)}. 


This is in EQP, by Lemma 3.1(5). 
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